Single-cell profiling of MC1R-inhibited melanocytes.

The human red hair color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to define a MC1R-inhibited Gene Signature (MiGS) comprising a large set of previously unidentified genes which may be implicated in melanogenesis and oncogenic transformation. We show that one of the candidate MiGS genes, TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into further mechanisms by which melanocytes with reduced MC1R signaling may regulate pigmentation and offer new candidates of study toward understanding how individuals with the RHC phenotype are predisposed to melanoma.

The MFSD12 p.Tyr182His common variant is sufficient to alter mouse agouti coat color.

MFSD12 functions as a transmembrane protein required for import of cysteine into melanosomes and lysosomes. The MFSD12 locus has been associated with phenotypic variation in skin color across African, Latin American, and East Asian populations. The frequency of a particular MFSD12 coding variant, rs2240751 (MAF = 0.08), has been reported to correlate with solar radiation and occur at highest frequency in Peruvian (PEL MAF = 0.48) and Han Chinese (CHB MAF = 0.40) populations, suggesting it could be causative for associated phenotypic variation in skin color. We have generated a mouse knock-in allele, Mfsd12Y182H , to model the human missense p.Tyr182His human variant. We demonstrate that the variant transcript is stably expressed and that agouti mice homozygote for the variant allele are viable with an altered coat color. This in vivo data confirms that the MFSD12 p.Tyr182His variant functions as a hypomorphic allele sufficient to alter mammalian pigmentation.

Haplotype-based analysis resolves missing heritability in oculocutaneous albinism type 1B.

Oculocutaneous albinism (OCA) is a rare disorder of pigment production. Affected individuals have variably decreased global pigmentation and visual-developmental changes that lead to low vision. OCA is notable for significant missing heritability, particularly among individuals with residual pigmentation. Tyrosinase (TYR) is the rate-limiting enzyme in melanin pigment biosynthesis and mutations that decrease enzyme function are one of the most common causes of OCA. We present the analysis of high-depth short-read TYR sequencing data for a cohort of 352 OCA probands, ∼50% of whom were previously sequenced without yielding a definitive diagnostic result. Our analysis identified 66 TYR single-nucleotide variants (SNVs) and small insertion/deletions (indels), 3 structural variants, and a rare haplotype comprised of two common frequency variants (p.Ser192Tyr and p.Arg402Gln) in cis-orientation, present in 149/352 OCA probands. We further describe a detailed analysis of the disease-causing haplotype, p.[Ser192Tyr; Arg402Gln] ("cis-YQ"). Haplotype analysis suggests that the cis-YQ allele arose by recombination and that multiple cis-YQ haplotypes are segregating in OCA-affected individuals and control populations. The cis-YQ allele is the most common disease-causing allele in our cohort, representing 19.1% (57/298) of TYR pathogenic alleles in individuals with type 1 (TYR-associated) OCA. Finally, among the 66 TYR variants, we found several additional alleles defined by a cis-oriented combination of minor, potentially hypomorph-producing alleles at common variant sites plus a second, rare pathogenic variant. Together, these results suggest that identification of phased variants for the full TYR locus are required for an exhaustive assessment for potentially disease-causing alleles.

Loss of MC1R signaling implicates TBX3 in pheomelanogenesis and melanoma predisposition.

The human Red Hair Color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA-sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to reveal a Pheomelanin Gene Signature (PGS) comprising genes implicated in melanogenesis and oncogenic transformation. We show that TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, is part of the PGS and binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into mechanisms by which MC1R signaling regulates pigmentation and how individuals with the RHC phenotype are predisposed to melanoma.

Correlation of age of onset and clinical severity in Niemann-Pick disease type C1 with lysosomal abnormalities and gene expression.

Niemann-Pick disease type C1 (NPC1) is a rare, prematurely fatal lysosomal storage disorder which exhibits highly variable severity and disease progression as well as a wide-ranging age of onset, from perinatal stages to adulthood. This heterogeneity has made it difficult to obtain prompt diagnosis and to predict disease course. In addition, small NPC1 patient sample sizes have been a limiting factor in acquiring genome-wide transcriptome data. In this study, primary fibroblasts from an extensive cohort of 41 NPC1 patients were used to validate our previous findings that the lysosomal quantitative probe LysoTracker can be used as a predictor for age of onset and disease severity. We also examined the correlation between these clinical parameters and RNA expression data from primary fibroblasts and identified a set of genes that were significantly associated with lysosomal defects or age of onset, in particular neurological symptom onset. Hierarchical clustering showed that these genes exhibited distinct expression patterns among patient subgroups. This study is the first to collect transcriptomic data on such a large scale in correlation with clinical and cellular phenotypes, providing a rich genomic resource to address NPC1 clinical heterogeneity and discover potential biomarkers, disease modifiers, or therapeutic targets.

Transcriptome of HPßCD-treated Niemann-Pick disease type C1 cells highlights GPNMB as a biomarker for therapeutics.

The rare, fatal neurodegenerative disorder Niemann-Pick disease type C1 (NPC1) arises from lysosomal accumulation of unesterified cholesterol and glycosphingolipids. These subcellular pathologies lead to phenotypes of hepatosplenomegaly, neurological degeneration and premature death. The timing and severity of NPC1 clinical presentation is extremely heterogeneous. This study analyzed RNA-Seq data from 42 NPC1 patient-derived, primary fibroblast cell lines to determine transcriptional changes induced by treatment with 2-hydroxypropyl-ß-cyclodextrin (HPßCD), a compound currently under investigation in clinical trials. A total of 485 HPßCD-responsive genes were identified. Pathway enrichment analysis of these genes showed significant involvement in cholesterol and lipid biosynthesis. Furthermore, immunohistochemistry of the cerebellum as well as measurements of plasma from Npc1m1N null mice treated with HPßCD and adeno-associated virus gene therapy suggests that one of the identified genes, GPNMB, may serve as a useful biomarker of treatment response in NPC1 disease. Overall, this large NPC1 patient-derived dataset provides a comprehensive foundation for understanding the genomic response to HPßCD treatment.

A custom capture sequence approach for oculocutaneous albinism identifies structural variant alleles at the OCA2 locus.

Oculocutaneous albinism (OCA) is a heritable disorder of pigment production that manifests as hypopigmentation and altered eye development. Exon sequencing of known OCA genes is unsuccessful in producing a complete molecular diagnosis for a significant number of affected individuals. We sequenced the DNA of individuals with OCA using short-read custom capture sequencing that targeted coding, intronic, and noncoding regulatory regions of known OCA genes, and genome-wide association study-associated pigmentation loci. We identified an OCA2 complex structural variant (CxSV), defined by a 143 kb inverted segment reintroduced in intron 1, upstream of the native location. The corresponding CxSV junctions were observed in 11/390 probands screened. The 143 kb CxSV presents in one family as a copy number variant duplication for the 143 kb region. In the remaining 10/11 families, the 143 kb CxSV acquired an additional 184 kb deletion across the same region, restoring exons 3-19 of OCA2 to a copy-number neutral state. Allele-associated haplotype analysis found rare SNVs rs374519281 and rs139696407 are linked with the 143 kb CxSV in both OCA2 alleles. For individuals in which customary molecular evaluation does not reveal a biallelic OCA diagnosis, we recommend preliminary screening for these haplotype-associated rare variants, followed by junction-specific validation for the OCA2 143 kb CxSV.

Genetic background modifies phenotypic severity and longevity in a mouse model of Niemann-Pick disease type C1.

Niemann-Pick disease type C1 (NPC1) is a rare, fatal neurodegenerative disorder characterized by lysosomal accumulation of unesterified cholesterol and glycosphingolipids. These subcellular pathologies lead to phenotypes of hepatosplenomegaly, neurological degeneration and premature death. NPC1 is extremely heterogeneous in the timing of clinical presentation and is associated with a wide spectrum of causative NPC1 mutations. To study the genetic architecture of NPC1, we have generated a new NPC1 mouse model, Npc1em1PavNpc1em1Pav/em1Pav mutants showed notably reduced NPC1 protein compared to controls and displayed the pathological and biochemical hallmarks of NPC1. Interestingly, Npc1em1Pav/em1Pav mutants on a C57BL/6J genetic background showed more severe visceral pathology and a significantly shorter lifespan compared to Npc1em1Pav/em1Pav mutants on a BALB/cJ background, suggesting that strain-specific modifiers contribute to disease severity and survival. QTL analysis for lifespan of 202 backcross N2 mutants on a mixed C57BL/6J and BALB/cJ background detected significant linkage to markers on chromosomes 1 and 7. The discovery of these modifier regions demonstrates that mouse models are powerful tools for analyzing the genetics underlying rare human diseases, which can be used to improve understanding of the variability in NPC1 phenotypes and advance options for patient diagnosis and therapy.This article has an associated First Person interview with the first author of the paper.

NPC1 Deficiency in Mice is Associated with Fetal Growth Restriction, Neonatal Lethality and Abnormal Lung Pathology.

The rare lysosomal storage disorder Niemann-Pick disease type C1 (NPC1) arises from mutation of NPC1, which encodes a lysosomal transmembrane protein essential for normal transport and trafficking of cholesterol and sphingolipids. NPC1 is highly heterogeneous in both clinical phenotypes and age of onset. Previous studies have reported sub-Mendelian survival rates for mice homozygous for various Npc1 mutant alleles but have not studied the potential mechanisms underlying this phenotype. We performed the first developmental analysis of a Npc1 mouse model, Npc1em1Pav, and discovered significant fetal growth restriction in homozygous mutants beginning at E16.5. Npc1em1Pav/em1Pav mice also exhibited cyanosis, increased respiratory effort, and over 50% lethality at birth. Analysis of neonatal lung tissues revealed lipid accumulation, notable abnormalities in surfactant, and enlarged alveolar macrophages, suggesting that lung abnormalities may be associated with neonatal lethality in Npc1em1Pav/em1Pav mice. The phenotypic severity of the Npc1em1Pav model facilitated this first analysis of perinatal lethality and lung pathology in an NPC1 model organism, and this model may serve as a useful resource for developing treatments for respiratory complications seen in NPC1 patients.

Identification of Gene Variants Associated with Melanocyte Stem Cell Differentiation in Mice Predisposed for Hair Graying.

Age-related hair graying is caused by malfunction in the regenerative potential of the adult pigmentation system. The retention of hair color over the life of an organism is dependent on the ability of the melanocyte stem cells and their progeny to produce pigment each time a new hair grows. Age-related hair graying is variable in association with genetic background suggesting that quantitative trait loci influencing this trait can be identified. Identification of these quantitative trait loci may lead to the discovery of novel and interesting genes involved in stem cell biology and/or melanogenesis. With this in mind we developed previously a sensitized, mouse modifier screen and discovered that the DBA/1J background is particularly resistant to melanocyte stem cell differentiation in comparison to the C57BL/6J background. Melanocyte stem cell differentiation generally precedes hair graying and is observed in melanocyte stem cells with age. Using quantitative trait loci analysis, we have now identified three quantitative trait loci on mouse chromosomes 7, 13, and X that are associated with DBA/1J-mediated variability in melanocyte stem cell differentiation. Taking advantage of publicly-available mouse sequence and variant data, in silico protein prediction programs, and whole genome gene expression results we describe a short list of potential candidate genes that we anticipate to be involved in melanocyte stem cell biology in mice.

A curated gene list for expanding the horizons of pigmentation biology.

Over the past century, studies of human pigmentary disorders along with mouse and zebrafish models have shed light on the many cellular functions associated with visible pigment phenotypes. This has led to numerous genes annotated with the ontology term "pigmentation" in independent human, mouse, and zebrafish databases. Comparisons among these datasets revealed that each is individually incomplete in documenting all genes involved in integument-based pigmentation phenotypes. Additionally, each database contained inherent species-specific biases in data annotation, and the term "pigmentation" did not solely reflect integument pigmentation phenotypes. This review presents a comprehensive, cross-species list of 650 genes involved in pigmentation phenotypes that was compiled with extensive manual curation of genes annotated in OMIM, MGI, ZFIN, and GO. The resulting cross-species list of genes both intrinsic and extrinsic to integument pigment cells provides a valuable tool that can be used to expand our knowledge of complex, pigmentation-associated pathways.

In Niemann-Pick C1 mouse models, glial-only expression of the normal gene extends survival much further than do changes in genetic background or treatment with hydroxypropyl-beta-cyclodextrin.

The Npc1nmf164 allele of Npc1 provides a mouse model for Niemann-Pick disease type C1 (NPC1), a genetic disease known to have a widely variable phenotype. The transfer of the Npc1nmf164 mutation from the C57BL/6J inbred strain to the BALB/cJ inbred strain increased the mean lifespan from 117.8days to 153.1days, confirming that the severity of the NPC1 phenotype is strongly influenced by genetic background. The transfer of another Npc1 allele, Npc1nih, to this background also extended survival of the homozygotes indicating that the modifying effect of BALB/cJ is not limited to a single allele of Npc1. The increased longevity due to the BALB/cJ background did not map to a previously mapped modifier on chromosome 19, indicating the presence of additional genes impacting disease severity. The previously studied Glial Fibrillary Acidic Protein promoter-Npc1 cDNA transgene (GFAP-Npc1) which only expresses NPC1 in astrocytes further extended the lifespan of Npc1nmf164 homozygotes on a BALB/cJ background (up to 600days). Hydroxypropyl-ß-cyclodextrin (HPßCD) treatment, not previously tested in the Npc1nmf164 mutant, extended life in the Npc1nmf164 homozygotes but not the transgenic, Npc1nmf164 mice on the BALB/cJ background. In all cases, lack of weight gain and early cerebellar symptoms of loss of motor control were found. At termination, the one mouse sacrificed for histological studies showed severe, diffuse pulmonary alveolar proteinosis suggesting that pulmonary abnormalities in NPC1 mouse models are not unique to the Npc1nih allele.

Loci associated with skin pigmentation identified in African populations.

Despite the wide range of skin pigmentation in humans, little is known about its genetic basis in global populations. Examining ethnically diverse African genomes, we identify variants in or near SLC24A5, MFSD12, DDB1, TMEM138, OCA2, and HERC2 that are significantly associated with skin pigmentation. Genetic evidence indicates that the light pigmentation variant at SLC24A5 was introduced into East Africa by gene flow from non-Africans. At all other loci, variants associated with dark pigmentation in Africans are identical by descent in South Asian and Australo-Melanesian populations. Functional analyses indicate that MFSD12 encodes a lysosomal protein that affects melanogenesis in zebrafish and mice, and that mutations in melanocyte-specific regulatory regions near DDB1/TMEM138 correlate with expression of ultraviolet response genes under selection in Eurasians.

BRG1 interacts with SOX10 to establish the melanocyte lineage and to promote differentiation.

Mutations in SOX10 cause neurocristopathies which display varying degrees of hypopigmentation. Using a sensitized mutagenesis screen, we identified Smarca4 as a modifier gene that exacerbates the phenotypic severity of Sox10 haplo-insufficient mice. Conditional deletion of Smarca4 in SOX10 expressing cells resulted in reduced numbers of cranial and ventral trunk melanoblasts. To define the requirement for the Smarca4 -encoded BRG1 subunit of the SWI/SNF chromatin remodeling complex, we employed in vitro models of melanocyte differentiation in which induction of melanocyte-specific gene expression is closely linked to chromatin alterations. We found that BRG1 was required for expression of Dct, Tyrp1 and Tyr, genes that are regulated by SOX10 and MITF and for chromatin remodeling at distal and proximal regulatory sites. SOX10 was found to physically interact with BRG1 in differentiating melanocytes and binding of SOX10 to the Tyrp1 distal enhancer temporally coincided with recruitment of BRG1. Our data show that SOX10 cooperates with MITF to facilitate BRG1 binding to distal enhancers of melanocyte-specific genes. Thus, BRG1 is a SOX10 co-activator, required to establish the melanocyte lineage and promote expression of genes important for melanocyte function.

Highly Efficient Cpf1-Mediated Gene Targeting in Mice Following High Concentration Pronuclear Injection.

Cpf1 has emerged as an alternative to the Cas9 RNA-guided nuclease. Here we show that gene targeting rates in mice using Cpf1 can meet, or even surpass, Cas9 targeting rates (approaching 100% targeting), but require higher concentrations of mRNA and guide. We also demonstrate that coinjecting two guides with close targeting sites can result in synergistic genomic cutting, even if one of the guides has minimal cutting activity.

Genomic analysis reveals distinct mechanisms and functional classes of SOX10-regulated genes in melanocytes.

SOX10 is required for melanocyte development and maintenance, and has been linked to melanoma initiation and progression. However, the molecular mechanisms by which SOX10 guides the appropriate gene expression programs necessary to promote the melanocyte lineage are not fully understood. Here we employ genetic and epigenomic analysis approaches to uncover novel genomic targets and previously unappreciated molecular roles of SOX10 in melanocytes. Through global analysis of SOX10-binding sites and epigenetic characteristics of chromatin states, we uncover an extensive catalog of SOX10 targets genome-wide. Our findings reveal that SOX10 predominantly engages 'open' chromatin regions and binds to distal regulatory elements, including novel and previously known melanocyte enhancers. Integrated chromatin occupancy and transcriptome analysis suggest a role for SOX10 in both transcriptional activation and repression to regulate functionally distinct classes of genes. We demonstrate that distinct epigenetic signatures and cis-regulatory sequence motifs predicted to bind putative co-regulatory transcription factors define SOX10-activated and SOX10-repressed target genes. Collectively, these findings uncover a central role of SOX10 as a global regulator of gene expression in the melanocyte lineage by targeting diverse regulatory pathways.

The transcription factors Ets1 and Sox10 interact during murine melanocyte development.

Melanocytes, the pigment-producing cells, arise from multipotent neural crest (NC) cells during embryogenesis. Many genes required for melanocyte development were identified using mouse pigmentation mutants. The variable spotting mouse pigmentation mutant arose spontaneously at the Jackson Laboratory. We identified a G-to-A nucleotide transition in exon 3 of the Ets1 gene in variable spotting, which results in a missense G102E mutation. Homozygous variable spotting mice exhibit sporadic white spotting. Similarly, mice carrying a targeted deletion of Ets1 exhibit hypopigmentation; nevertheless, the function of Ets1 in melanocyte development is unknown. The transcription factor Ets1 is widely expressed in developing organs and tissues, including the NC. In the chick, Ets1 is required for the expression of Sox10, a transcription factor critical for the development of various NC derivatives, including melanocytes. We show that Ets1 is required early for murine NC cell and melanocyte precursor survival in vivo. Given the importance of Ets1 for Sox10 expression in the chick, we investigated a potential genetic interaction between these genes by comparing the hypopigmentation phenotypes of single and double heterozygous mice. The incidence of hypopigmentation in double heterozygotes was significantly greater than in single heterozygotes. The area of hypopigmentation in double heterozygotes was significantly larger than would be expected from the addition of the areas of hypopigmentation of single heterozygotes, suggesting that Ets1 and Sox10 interact synergistically in melanocyte development. Since Sox10 is also essential for enteric ganglia development, we examined the distal colons of Ets1 null mutants and found a significant decrease in enteric innervation, which was exacerbated by Sox10 heterozygosity. At the molecular level, Ets1 was found to activate an enhancer critical for Sox10 expression in NC-derived structures. Furthermore, enhancer activation was significantly inhibited by the variable spotting mutation. Together, these results suggest that Ets1 and Sox10 interact to promote proper melanocyte and enteric ganglia development from the NC.

Ectopic differentiation of melanocyte stem cells is influenced by genetic background.

Hair graying in mouse is attributed to the loss of melanocyte stem cell function and the progressive depletion of the follicular melanocyte population. Single-gene, hair graying mouse models have pointed to a number of critical pathways involved in melanocyte stem cell biology; however, the broad range of phenotypic variation observed in human hair graying suggests that additional genetic variants involved in this process may yet be discovered. Using a sensitized approach, we ask here whether natural genetic variation influences a predominant cellular mechanism of hair graying in mouse, melanocyte stem cell differentiation. We developed an innovative method to quantify melanocyte stem cell differentiation by measuring ectopically pigmented melanocyte stem cells in response to the melanocyte-specific transgene Tg(Dct-Sox10). We make the novel observation that the production of ectopically pigmented melanocyte stem cells varies considerably across strains. The success of sensitizing for melanocyte stem cell differentiation by Tg(Dct-Sox10) sets the stage for future investigations into the genetic basis of strain-specific contributions to melanocyte stem cell biology.

Method of euthanasia influences the oocyte fertilization rate with fresh mouse sperm.

In vitro fertilization (IVF) is used to produce mouse embryos for a variety of reasons. We evaluated the effect of the method of euthanasia on the fertilization rate in 2 different IVF protocols. Oocytes collected from C57BL/6J female mice euthanized by CO2 inhalation or cervical dislocation were used in IVF with fresh sperm from either wild-type or genetically engineered C57BL/6J. Compared with CO2 inhalation, cervical dislocation improved the resulting rate of fertilization by 18% in an IVF method using Cook media and by 13% in an IVF method using methyl-B cyclodextrin and reduced glutathione. The lower fertilization rate due to euthanasia by CO2 inhalation was accompanied by changes in blood pH and body temperature despite efforts to minimize temperature drops. In our hands, euthanasia by cervical dislocation improved fertilization rates and consequently reduced the number of egg-donor mice required.